Atomic ordering effects and bond indices in quaternary systems Si-(Be,Al)-0-N 
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We propose a microscopic method for modeling atomic ordering effects in multi-component solid 
solutions, which is based on minimizing the number of unfavorable interatomic bonds in the system. 
Atomic ordering effects are studied in the solid solutions Sie-i^Alj^OajNg-j, (/3-SiAlONs) and hypo- 
thetical Si6-a;Bei,02i:N8-2a; iu the /3-Si3N4-Be-0 system. It is established that in /3-SiAlONs, A1,0 
atoms form separate quasi-one-dimensional "channels". On the contrary, in Si6-i^Bea;02iN8-2i:, Be 
and O atoms constitute a whole cluster. This means that homogeneous solid solutions cannot be 
formed in the /3-Si3N4-Be-0 system under equilibrium conditions. 
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I. INTRODUCTION 

Atomic ordering effects (AOE) have been discovered in 
a wide range of non-stoichiometric and multi-component 
alloys, solid solutions (SS) and compounds. The influence 
they produce on the properties is great and comparable 
with the effect of altering the chemical composition of 
materials 

In the majority of theoretical simulations of multi- 
component materials, atomic ordering effects are studied 
using integral energy characteristics (for example, free, 
total or cohesion energies) of a limited number of atomic 
configurations of a system set a priory [^j. Generally, 
this approach is applied to binary or ternary compounds, 
whereas for more complex systems its application is con- 
fined by computer power due to rapidly increasing sizes 
of cells to be considered and the number of atomic con- 
figurations. 

In the present report, we describe a new method for 
studying AOE based on comparative analysis of detached 
types of interatomic bonds in multi-component systems. 
In the framework of this method, optimal atomic configu- 
rations are chosen by minimizing unfavorable interatomic 
bonds. 

The method is illustrated using quaternary solid solu- 
tions Si6-sAla;02,N8_x (/3-SiAlONs) considered earlier in 
ref. 1^ as an example. It is also applied to predict atomic 
ordering effects in the hypothetical /3-Si3N4-based solid 
solutions Si6-iBej;02xN8-2x- 



II. THE METHOD FOR MODELING ATOMIC 
ORDERING EFFECTS: /3-SIALONS 

As is known, the stability of a thermodynamic system 
corresponds to the free energy or, at zero temperature, 
total energy minimum condition. In terms of the chem- 
ical bonding theory, this condition implies the optimal 



type of bands filling, when all bonding states are occu- 
pied and all antibonding states are vacant. 

In the proposed method, AOE are determined by 
searching for the multi-component system configura- 
tions having the maximum number of the most favor- 
able interatomic bonds. To characterize the detached 
type of paired bonds, crystal orbital overlap populations 
(COOP) and the numbers of filled antibonding states 
(FAS) may be used. In this work, these values are cal- 
culated employing the MuUiken analysis by the tight- 
binding band method with matrix elements parameter- 
ization according to the Hiickel theory Other tech- 
niques of bond characterization can be used as well. 

The application of this approach can be demonstrated 
on the four-component SS in the Si-Al-O-N system with 
the basic structure /3-Si3N4 (the so-called /?-SiA10Ns of 
the formal stoichiometry Sie-ajA^OajNg-j,, where x = 
0-4.2). 

A previous investigation of /3-SiAlONs by the tight- 
binding band method showed that Al and O atoms 
form extended quasi-one-dimensional (ID) structures 
constituted by 12-atomic rings — the so-called "impu- 
rity channels", fig. |l|. Subsequent studies by the DFT 
cluster [ H and the band structure full-potential LMTO 
methods showed that ordered structures exhibited the 
best stability. Besides, their chemical bonding and elec- 
tronic properties were examined. 

A detailed analysis of bond indices in /3-SiAlONs 
demonstrated Q that the most unfavorable bonds are Si- 
O, which contain the maximum number of FAS. The an- 
tibonding states of Al-0 bonds are less populated. Si-N 
and Al-N bonds do not contain any FAS, while COOP(Si- 
N= 0.58 e)>COOP(Al-N= 0.24 e). Hence, a bonds hi- 
erarchy taking into account their stabilizing effects can 
be represented by the so-called sequence of bonds prefer- 
ences: 



Psi_o < -Pai-o < ^Al-N < Pi 



Si-N, 



(1) 
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where Pa-b refiects the preference of A-B-type bonds 
presence in the system. 

Thus, in order to determine the superstructure of a 
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FIG. 1: The impurity channels in the structure of 
Sis-xALO^Ns-x (x = 1.3). 

multi-component material, it is necessary to choose such 
distributions of atoms over lattice sites, in which the 
number of bonds is minimized consecutively in accor- 
dance with (|l]). 

For example, for /3-SiAlONs at the first step the super- 
structures with the minimum number of Si-0 bonds are 
chosen from all possible ones according to (j^). Among 
them we choose the superstructures with the minimum 
number of Al-0 bonds etc. As a result, in a system of a 
given composition, the superstructures are found, which 
possess the same number of the same types of bonds, but 
differ only by their positions in the supercell. 

To determine the atomic order type, the following al- 
gorithm is applied. The SS Sie-xAljjOajNg-a; of a cer- 
tain composition modelled by a supercell is considered. 
In this work we make use of the 126-atomic supercells 
(3x3x1) analogous to those applied in ||^. For example, 
the composition of the SS with x = 2.0 (Si4Al202N6) 
is described by the Si36Ali80i8N54 cell. The following 
items are fulfilled: 

(i) Atoms (Si,Al,N,0) are randomly distributed over the 

supercell sites so that the result does not depend 
on initial conditions. 

(ii) A random substitution, which is the only permitted 

one for this system (for example Si-^Al, N^O), 
is performed, while all other types of substitution 
(Si— >0, N— >A1) are not considered. 

(iii) If the correspondence to (Q) improves (the number 
of favorable bonds increases and the number of un- 
favorable bonds3 decreases), the step is accepted, 
otherwise it is rejected. 

(iv) Items (ii) and (iii) are repeated until the number of 
steps denied in succession does not exceed iVi. 

(v) Items (i-iv) are repeated until the number of struc- 

tures denied in succession does not exceed N2 ■ This 
cycle is necessary since the result of cycle (iv) de- 
pends on the initial position of atoms and on the 



first accepted substitution (item iii) defining the 
subsequent shape of the superstructure. 

The choice of A^i, N2 values is determined by the su- 
percell size and composition. In the present work they 
range from 10000 to 100000 and from 10 to 400 respec- 
tively. 

An application of the above algorithm for /3-SiAlONs 
with the same cells and compositions as in Q shows that 
the resulting superstructures shapes (ID impurity chan- 
nels, fig. |l|) are in complete agreement with those ob- 
tained in a series of complicated calculations of total en- 
ergies. The formation of impurity channels in /3-SiAlONs 
leads to the disappearance of the most unfavorable Si-0 
bonds, while the number of Al-0 bonds is minimal. 

The evident drawback of the aforesaid scheme con- 
cerns condition (|l|) , which reflects only the effect of paired 
bonds A-B. It should not be excluded that in some cases 
the formation of A„B„i clusters is more preferable in spite 
of low preference of A-B bonds. This restriction may be 
removed by using a sequence of clusters preferences in- 
stead of a sequence of bonds preferences. In respect to 
/3-SiAlONs, the above restriction does not permit differ- 
entiating between distant and adjacent positions of im- 
purity channels. 



III. ELECTRONIC STRUCTURE AND ATOMIC 
ORDERING EFFECTS IN THE /5-SI3N4-BE-0 
SYSTEM 

A. Bond indices and chemical composition 

At the first stage of modeling the hypothetical SS in 
the quaternary system /3-Si3N4-Be-0, we analyzed the 
electronic properties and bond indices of pure 5-Si3N4 
and that containing Be and O impurities, fig. ^. It is 
seen that in the binary nitride /3-Si3N4 (cell Si54N72, va- 
lence electrons concentration (VEC) equals 576 e/cell), 
all antibonding states are vacant, while all bonding states 
are completely occupied. 

O^N or Be— >Si substitutions (cells Si540N7i and 
Si53BeN72 respectively) reduce the stability of the (3- 
Si3N4:0 and /?-Si3N4:Be systems. In the first case this is 
due to partial filling of antibonding states when VEC in- 
creases (577 e/cell), while in the second case — due to a 
depopulation some bonding states when VEC decreases 
(574 e/cell). A "mutual compensation" of the destabi- 
lizing effects of single impurities Be, O can be achieved 
by choosing a composition, at which the VEC remains 
576 e/cell, i.e. when the impurities are introduced in 
the proportion (Be-|-20). In this case the type of bands 
filling (as compared with the basic /3-Si3N4) does not 
change and the number of filled antibonding and empty 
bonding states is minimal, fig. ^j. Thus, the formal 
stoichiometry of the SS in the /3-Si3N4-Be-0 system is 

Si6-i;Bea;022;N8-2a:- 
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FIG. 2: Crystal orbital overlap populations of bonds in Si3N4 (a) and impurity systems Si3N4:0 (b), Si3N4:Be (c), 
Si3N4:(Be+20) (d) simulated by Si54N72, Si540N7i, Si53BeN72 Si53Be02N7o cells, respectively. 



B. Short-range atomic ordering 

When heterovalent substitutions occur in /3-Si3N4, the 
adjacent arrangement of donour (O) and acceptor (Be) 
atoms (i.e. the formation of {Be02} complexes) is more 
probable, otherwise the charge screening effect hinders 
the above mentioned compensation of bonds fiUing. 



C. Long-range atomic ordering 

Analysis of separate bonds in the Si6-a:Bea;02a;N8_2a; 
system shovifs that Si-0 and Be-0 bonds consist of filled 
antibonding states, and their number in Si-0 is greater 
than in Be-0 (fig. Determining atomic ordering ef- 
fects by the proposed algorithm, we search for atomic 
configurations satisfying the condition of type (||): 



(2) 



In other words, for the considered series of solid solu- 
tions of variable concentrations Si6-a;Bea;02£cN8_2a; {x ~ 
— 0.5, supercells Si53Be02N7o - Si24Be3o06oNi2), we 
determine the configurations of modelled cells, in which 
the number of less advantageous bonds Si-0 and Be-0 is 
successively minimized. 

The results of the simulation are depicted in fig. ^ and 
table |. Their major difference from pSie-x-^^xOx^s-x 
is the presence of Si-0 bonds, the number of which grows 
monotonously with the concentration of Be and O. This 
indicates that Si6-i^Be2:02a;N8-2a; is less stable in com- 
parison with Sig-xAUOxNg-a;- When x is small (cells 



Si53Be02N7o - Si5iBe306N66), the system contains Be- 
N bonds, and the number of "unfavorable" Si-0 bonds 
grows with x most rapidly. The noted peculiarities are 
due to a small number of Be and O atoms in the supercell 
and its limited size. In real crystals this situation is not 
possible. 

The smallest cluster, which does not contain any Be-N 
bonds, is made up of four {BCO2} complexes (supercell 
Si5oBe408N64, fig. ^). As the number of Be, O atoms 
grows further, the number of Si-0 bonds increases lin- 
early and equales the double number of {Be02} com- 
plexes in the cell. This fact shows that the system has 
no definite order because, for example, one cluster of 8 
complexes {Be02} contains the same set of bonds as two 



TABLE I: The numbers of different-type bonds in hypothet- 
ical SS Si6-i:Bea;02a;N8-2i: as a function of the modelled cell 
composition. 



Cell 


Bonds 


Si-0 


Be-0 


Be-N 


Si-N 


Si53BeO2N70 


3 


3 


1 


209 


Si52Be204N68 


5 


7 


1 


203 


Si5iBe306N66 


7 


11 


1 


197 


Si5oBe408N64 


8 


16 





192 


Si49Be50ioN62 


10 


20 





186 


Si48Be60i2N6o 


12 


24 





180 


Si47Be70l4N58 


14 


28 





174 


Si46Be80l6N56 


16 


32 





168 


Si26Be28 056Ni6 


56 


112 





48 


Si25Be29 058Nl4 


58 


116 





42 
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60 


120 





36 
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FIG. 3: Superstructures in cells of Si5oBe408N64 (a), Si46Be80i6N58 (b), Si43BeiiO22N50 (c,d), Si4iBei3 026 N46 (e), 
Si4oBei4028N44 (f) compositions. 



separate clusters each containing 4 {Be02} complexes 
(fig. H). As a result, atomic configurations at a definite 
supercell composition may have the form of extended, 
compact or multiple clusters constituted by {Be02} com- 
plexes, fig. ||. 

The analysis of bonds allows us to conclude that the 
formation of a whole compact cluster with the maximal 
size is most probable. In this case, the "depopulation" 
of FAS occurs analogously as in precipitated phases of 
concentration polytypes in the AlN-Si-0 system |j, 
where antibonding states are present only in "interlayer" 
bonds Al-0 and are absent within the "block of defects" . 

In Si6-xBea;022;N8_22:, an increase in the size of the 



"impurity cluster" brings about a growth of the num- 
ber of "intracluster" Be-0 bonds reducing the number of 
FAS. For example, COOP spectra of "intracluster" and 
"frontier" Be-0 bonds of the Si46Be80i6N56 supercell are 
depicted in fig. ^. It can be seen that intracluster bonds 
do not possess any antibonding states. This indicates 
that the formation of a compact cluster of the largest 
size is more advantageous. The growth of a cluster in 
Si6-xBex022;N8_2x is not limited meaning the formation 
of a heterophase state of the system (/3-Si3N4/BeO). 

Thus, the proposed semiempirical approach based on 
quantum-chemical analysis of separate interatomic bonds 
makes it possible to determine the most probable ef- 
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FIG. 4; COOP spectra of "intracluster" (a) and "frontier" (b) 
Be-O bonds in a cluster of Be, O impurities in Si43Beii022N5o 
(fig- |c). 

fects of short-range (clusterization) and long-range (su- 
perstructures) atomic ordering in multi-component solid 
solutions and to find out the factors responsible for their 
stability. 

Analysis of the differences between Si6_a;Be2;022;Ng_2a; 
and Si6_a;Ala,02;N8_a; solid solutions shows that the im- 
possibility of formation of homogeneous SS in the (3- 
Si3N4-Be-0 system is due to the excess of oxygen atoms 
with regard to those in /3-SiAlONs. This leads to the 
presence of Si-0 bonds in the system, which does not al- 
low the impurities to form compact structures of a limited 
size like impurity channels. Therefore it can be suggested 
that SS in the /3-Si3N4-Be-0 system may be stabilized 
by altering the stoichiometry of the compounds (compo- 
sition Sie-ccBejjOajNg-a;) or by partial replacement of Si 
atoms by elements of the IV-th group, for example, Zr 

(Si6_a;Bea;/2Zra:/20^N8_^). 

Undoubtedly, the detailed analysis of the potentially 



stable structures, which were suggested in the frame- 
work of the used approach, implies further investiga- 
tions of their energy states by ab-initio methods for the 
conformation of the presented results. At present, the 
above studies of Si(i^xQGx02x^%-2x, Sie-ajBe^jOajNg-a; 
and Si6-a;Bej./2Zrj./20a;N8_a; SS are being performed by 
the full-potential LMTO method. 



IV. CONCLUSIONS 

In this paper we proposed a new microscopic semi- 
empirical method for describing atomic ordering effects 
in multi-component systems. This method is based on 
the analysis of interatomic bonds of separate types and 
includes a procedure of minimizing the number of unfa- 
vorable bonds. 

The method is used to study AOE in solid solutions 
of variable compositions: Si6-2;Ala;02:N8_a; (/3-SiAlONs) 
and hypothetical Si6-a;Be2;02a;Ng_2a; in the /9-Si3N4-Be- 
O system. It is established that in /3-SiAlONs Al, O 
atoms form separate quasi-one-dimensional "impurity 
channels" . On the contrary, in Si6-a;Bea;02a;N8_2a; Be, 
O atoms trend to form a whole compact cluster. This 
implies that in the /3-Si3N4-Be-0 system homogeneous 
solid solutions cannot be formed under equilibrium con- 
ditions. 

Potential ways of stabilizing the SS in the /3-Si3N4- 
Be-0 system are supposed to be a variation of its sto- 
ichiometry (reduction in the oxygen concentration) or 
partial replacement of silicon by other atoms from the 
IV-th group, for example, zirconium. 
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